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ABSTRACT 


-Y 

There  are  three  atmospheric  processes  responsible  for  the  degrada¬ 
tion  of  tie  transmission  of  optical  images  and  electro-optical  energy: 
aerosol  extinction,  molecular  absorption,  and  turbulent  distortion  (scin¬ 
tillation  and  beam  wander).  As  a  part  of  the  Marine  Aerosol  Generation 
and  Transport  experiment  (MAGAT-80),  light  transmission  characteristics 
(refractive-index  structure  function  parameter,  Cjj2,  and  total  ex¬ 
tinction  coefficient,  a)  were  measured  optically  on  a  13.3  km  path  across 
Monterey  Bay.  C^2  and  a  can  also  be  calculated  from  micrometeoro- 
logical  data  (aerosol  spectra,  turbulence  and  mean  meteorological  para¬ 
meters).  This  report  is  a  compilation  of  the  preliminary  analysis  of 
path-averaged  (aircraft)  and  midpoint  (ship)  micrometeorological  data, 
including  calculations  of  the  relevant  optical  parameters  for  comparison 
with  the  optical  measurements. 
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A.  INTRODUCTION 


Light  propagating  through  the  atmosphere  is  not  only  scattered  and 
absorbed  by  aerosols  and  molecules,  but  the  wavefronts  are  deflected  and 
distorted  by  turbulence.  The  evaluation  and  application  of  optical, 
electro-optical,  and  laser  systems  requires  reliable  data  and  a  tested 
physical  model  of  these  atmospheric  effects. 

NFS  personnel  recently  participated  in  a  large-scale  field  experiment 
designed  to  iirprove  and  verify  certain  overwater  models  of  these  atmos¬ 
pheric  processes  for  the  U.S.  Navy.  The  experiment,  Marine  Aerosol 
Generation  and  Transport  (MAGAT) ,  was  the  brainchild  of  Professors 
Kenneth  L.  Davidson  and  Gordon  E.  Schacher  of  the  Environmental  Physics 
Group  at  the  Naval  Postgraduate  School  (NPS)  in  Monterey,  California. 
Other  government  installations  involved  in  the  planning  of  the  experiment 
were  the  Naval  Ocean  Systems  Center  and  the  Naval  Environmental  Predic¬ 
tion  Research  Facility.  The  Electro-Optics/Meteorology  (BO/MET)  Program, 
the  High  Qiergy  Laser  (HEL)  Program,  and  the  Naval  Air  Systems  Command 
provided  funding  for  the  project.  MAGAT  was  held  from  April  28  to  May  9, 
1980,  in  the  vicinity  of  Monterey  Bay. 

The  first  phase  of  the  experiment  dealt  with  the  ccrpatibility  of 
optical  and  micrometeorological  propagation  theory.  In  crop  .ration  with 
the  NPS  Physics  Department,  Optical  Propagation  Group,  direct  measure¬ 
ments  of  optical  extinction  and  scintillation  across  Monterey  Bay  were 
conpared  with  both  marine  surface  layer  model  predictions  and  aerosol  and 
turbulence  data  obtained  at  the  midpoint  of  the  optical  path  from  the 
Research  Vessel  R/V  ACANIA.  In  addition,  path  averages  of  aerosols  and 
turbulence  were  obtained  by  flying  an  instrumented  aircraft 
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the  entire  length  of  the  13.3  kilometer  optical  path  at  altitudes  varying 
from  3.5  to  20  meters  above  the  sea  surface.  The  aircraft  measurements 
were  made  in  cooperation  with  Airborne  Research  Associates  of  Boston, 
Massachusetts.  This  report  is  a  preliminary  analysis  of  the  aircraft  and 
ship  measurements  for  the  first  phase. 

The  second  phase  of  the  experiment  involved  an  ambitious  attempt  to 
extend  dynamic  models  of  the  evolving  marine  atmospheric  boundary  layer 
to  include  aerosol  and  turbulence  profiles.  This  phase  of  the  experi¬ 
ment,  conducted  in  a  region  30  to  50  nautical  miles  off  the  coast  of  Mon¬ 
terey,  required  periodic  monitoring  of  aerosol  and  micrcmeteorological 
variables  from  the  surface  to  5  kilometers.  These  duties  were  shared  by 

the  aircraft  and  the  R/V  ACANIA  (which  utilized  various  remote  sensing 
* 

techniques) .  The  analysis  of  the  second  phase  of  the  aircraft  measure¬ 
ments  will  be  covered  in  a  separate  report. 

B.  BACKGROUND 

1.  Optical  Parameters 

The  two  atmospheric  optical  properties  of  primary  interest  are 
total  extinction  and  refractive-index  structure  function  parameter, 
the  extinction  has  several  aonpcnents:  molecular  scattering  and 
absorption  (B  *  Bg  +  and  aerosol  scattering  and  absorption  (a  3  as  +  aa )  • 
Thus,  the  extinction  parameterizes  the  loss  of  light  energy  as  it  is 
scattered  out  of  the  bean  or  absorbed  by  the  molecular  and  particulate 
constituents  of  the  atmosphere.  The  distortion  and  tilt  of  image  wave 
fronts  by  atmospheric  turbulence'  is  parameterized  by  C^2. 

We  can  write  C^2  as  a  function  of  tenperature  (Op2)  and 
water  vapor  (Cq2)  turbulence  structure  function  parameters 

Cjj2  *  (79  x  lo”6  PA2)2(Cj,2  +  0.113  Qjg  +  3.2  x  10'3  CQ2)  (1) 


2 


MM 


where  P  is  the  pressure  in  nb,  T  the  absolute  temperature  and  C^q  the 
temperature-humidity  oospectral  structure  function  parameter.  C^2  can 
be  obtained  in  three  ways:  1)  optical  measurement,  2)  measurement  of 
G r2,  Giq,  and  Cq2,  and  3)  calculation  of  Gj2,  Gjq,  and  Cq2  from  bulk 
meteorological  data  (water  temperature,  air  temperature,  humidity  and 
wind  speed) . 

The  total  extinction  (a  +  8)  can  be  measured  optically  by  deter¬ 
mining  the  reduction  in  beam  intensity  over  sene  suitable  optical  path. 
The  separate  components  can  be  calculated  from  micrometeorological  data. 
The  molecular  extinction  can  be  obtained  from  the  LOWTRAN  model  developed 
by  the  Air  Force  Geophysics  laboratory  (Selby  et.  al,  1978).  The  aerosol 

extinction  can  be  calcualted  from  the  aerosol  spectral  density,  N(r) . 

(<*>  2 

r  =  I  rrr  E(n,X)  N(r)  dr  (2) 

where  r  is  the  particle  radius,  E(n,X)  the  total  scattering  efficiency  at 
wavelength,  X ,  and  refractive-index,  n. 


2.  Turbulence  Scaling  Parameters 

Since  the  details  of  surface  layer  scaling  are  covered  in  pre¬ 
vious  reports  (Fairall,  et.  al.,  1980)  this  discussion  will  be  limited  to 
a  few  basic  definitions.  Near  the  surface,  the  height  above  the  surface, 
Z,  can  be  normalized  by  the  Monin-Cbukhov  stabiltiy  length,  L.  We  can 
then  represent  the  mdcromebeorological  properties  in  terms  of  scaling 
parameters  and  dimensionless  functions  of  f  *  Z/L, 

Qj,2  -  T*2  Z"2/3  f  (5)  (3a) 

CQ2  =  Q*2  Z_2/3  Af  (5)  (3b) 

Sj  -  t.q.z-2/3  aV2  f(C)  ^  (3c) 
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where  T*  and  Q*  are  the  temperature  and  humidity  seeding  parameters,  f (£) 
is  a  dimensionless  function  (Wyngaard  et.  ad.,  1971),  r^  is  the  tem¬ 
perature-humidity  correlation  parameter  (about  0.8)-  and  A  is  a  constant 
(about  0.6). 

The  rate  of  dissipation  of  turbulent  kinetic  energy,  ? ,  can  be 
similarly  represented 


e 


u* 


KZ 


g(C) 


(4) 


where  u*  is  the  friction  velocity  and  K  is  Von  Kantian' s  constant  (0.35). 
The  scaling  length  is  given  by 


L 


2 

T  u»  _ 

gK  (T*  +0.61  TQ*/P> 


(5) 


where  g  is  the  acceleration  of  gravity  and  p  is  the  density  of  air. 

Note  that  the  scaling  parameters  are  related  to  the  surface 
fluxes  of  momentum  (t  ■  Reynolds  stress),  temperature  (Qo)  and  water 
vapor  (Mo) 


2 

(6a) 

T  =  PU* 

Qo  -  -U*T* 

(6b) 

Mq  =  -u*Q* 

(6c) 

3.  Bulk  Parameterization 

Although  the  scaling  parameters  can  be  determined  from  either 
direct  flux  measurements  or  from  measurements  of  Cj.2,  Cq2,  and  e , 
the  difficulty  of  these  measurements  has  led  to  the  development  of  a 
method  that  utilizes  bulk  meteorological  quantities  (wind  speed,  u,  tem¬ 
perature,  T,  and  water  vapor  density,  Q) .  In  this  case,  the  scaling 
parameter  for  X  (X  *  u,T,Q)  is  obtained  from  the  difference  in  X  from 
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the  sea  surface  (Xs)  to  some  reference  height  (usually  10m)  in  the 
atmosphere. 

*.  -  Cx  *10  -  V  '7> 

where  is  the  drag  coefficient  for  x  (typically,  c*  =  1.3  x  10-3 
over  the  ocean) .  Further  details  on  the  bulk  method  can  be  found  in 
Davidson  et.  al.  (1980). 

C.  INSTROMEMTOTICN 

1.  Aerosol 

The  aerosol  spectra  were  measured  with  optical  particle  coun¬ 
ters  made  by  Particle  Measurement  Systems  (PMS)  of  Boulder,  Colorado. 

The  R/V  ACANIA  used  the  standard  NPS  system  aonsisting  of  two  probes, 
the  classical  scattering  (CSAS)  and  the  active  scattering  (ASAS) , 
controlled  by  a  DAS-32  with  computer  interfacing.  This  system  measures 
aerosols  in  90  size  channels  from  0.09u  to  14. 0y  radius.  The  aircraft 
aerosol  data  were  obtained  using  a  PMS  model  ASSAP  on  loan  from  NGSC. 

This  system  has  60  size  channels  from  0.28u  to  14. Qy  radius. 

2.  Aircraft  Meteorology 

The  aircraft  micrometeorological  parameters  are  logged  on  a 
computer  controlled  (HP9835)  twenty  channel  data  acquisition  system. 

Bach  parameter  is  sanpled  every  2.5  seconds  with  a  two-scan  average 
stored  every  5  seconds.  The  data  is  stored  an  magnetic  tape  cassette 
with  a  four  hour  capacity.  A  brief  description  of  the  mi  crane teoro- 
logical  data  is  given  in  Table  1.  Further  details  on  aircraft  instrumen¬ 
tation  can  be  found  in  Fairall  (1979). 
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D.  ANALYSIS 


1.  Aerosol 

The  aerosol  analysis  techniques  for  the  ship  and  aircraft  are 
basically  the  same.  The  N(r)  spectrum  is  calculated  for  half-hour  aver¬ 
ages  on  the  ship  and  path  averages  for  the  aircraft  (about  2  minutes). 

The  spectrum  is  fit  in  LOG(N(r)),  LOG(r)  space  with  a  seventh  order  poly¬ 
nomial  for  0.09p<r<7p,with  a  linear  fit  for  r>7jj.  The  extinction  is  cal¬ 
culated  using  these  fits  for  0. 03y<r<30ii  on  the  ship  and  for  0.1y<r<15,j 
for  the  aircraft.  This  calculation  is  discussed  in  depth  in  Schacher  et. 
al.  (1980). 

The  method  was  developed  for  the  ship  system  and  adapted  for 
use  with  the  aircraft.  Because  of  the  greater  statistical  scatter  in  the 
N(r)  spectrum  from  the  aircraft  probe,  the  polynomial  fit  is  subject  to 
occasional  "instabilities".  Should  this  occur,  the  polynomial  fit  will 
bear  no  resemblance  to  the  N(r)  data.  Another  symptom  of  this  instabil¬ 
ity  is  the  occurrence  of  large  polynomial  coefficients.  Due  to  the  pre¬ 
liminary  nature  of  this  report,  the  data  have  been  left  unedited.  The 
reader  is  cautioned  to  use  common  sense  when  attempting  to  use  these  re¬ 
sults. 

The  aircraft  and  ship  aerosol  extinctions  were  compared  in  a 
series  of  flybys.  Since  the  ship  system  is  newer,  has  a  wider  range, 
better  sensitivity  and  is  better  understood,  we  decided  to  correct  the 
aircraft  extinctions  to  agree  with  the  ship.  The  correction  factors  are 
given  in  Table  II. 


TABLE  I.  Aircraft  Meteorological  Data 


Channel  Data  Symbol  Sensor 


1 

Pressure 

P 

National  Semiconductor 

2 

Temperature 

T 

Platinun  resistor 

3 

Temperature 

T 

Vortex  (NRL) 

4 

Dew  Point 

Td 

Cooled  mirror 

5 

Sea  Surface  T 

TS 

PRT-5  (IR) 

6 

Electric  Field 

E 

Radioactive  probe 

7 

- 

8 

Refractivity 

N 

Microwave  cavity  (NAC) 

9 

Water  Vapor  Density  Q 

Lyman-a,  mean  (NRL) 

10 

Air  Speed 

0 

Hot  wire,  mean 

11 

Dissipation 

Hot  wire,  fluctuation 

12 

N  structure  funct. 

%2 

Microwave,  fluctuation 

13 

T  structure  funct. 

Cr2 

Microthermal,  fluctuation 

14 

Q  structure  funct. 

Cq2 

Lyman-oc,  fluctuation 

15 

- 

16 

- 

17 

- 

18 

- 

19 

Electrical  Conduc¬ 
tivity 

X 

Flat  plate 

20 

— 
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TABLE  II.  Ratio  of  Ship  to  Aircraft 
Extinction  Coefficient  Values 


Wavelength ,  um 

Before  5/4/80 

After  5/3/80 

0.63 

3.8 

1.8 

0.84 

5.2 

2.0 

1.06 

7.1 

2.5 

These  factors  are  based  only  on  the  open  ocean  aortparisons .  The  Monterey 
Bay  comparisons  were  not  included  so  that  the  ship  and  aircraft  optical 
conparisons  could  be  considered  independent.  The  complete  set  of  correc¬ 
tion  factors  is  shown  in  Fig.  1. 

2.  Micrcmeteorology 

The  methods  and  equations  used  to  obtain  the  basic  parameters 
given  in  Table  I  have  been  described  in  Fairall  (1979).  Once  these 
meteorological  parameters  are  in  hand,  one  can  calculate  the  scaling 
parameters  (Section  B2)  using  either  turbulence  or  bulk  quantities.  Since 
we  did  not  have  mean  wind  speed  available  for  the  aircraft,  we  did  bulk 
calculations  using  a  hybrid  method  where  the  dissipation  rate,  e,  is  used  to 
obtain  u*  (Eq.  4) . 

3.  Optical  Data 

The  optical  extinction  ceff icients,  as  obtained  from  the  optics 
group,  represent  total  extinction  due  to  aerosols  and  air  molecules.  The 
molecular  components  were  calculated  using  DCWTRAN  II IB  and  subtracted 
from  the  total  to  leave  only  the  aerosol  extinction.  The  ICWTRAN  values 
used  are  given  below 

Wavelength,  um  0.63  0.84  1.03  1.06 

Molecular  8,  tan"1  0.01  0.04  0.00  0.00 

A  description  of  the  optical  measurements  is  given  by  Crittenden  et.  al. 
(1980). 
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Ratio  of  ship  to  aircraft  aerosol  extinction  coefficients.  The 
horizontal  lines  represent  the  factors  used  to  correct  the  aircraft 
values  to  agree  with  the  ship. 


E.  RESULTS 


The  measurements  were  made  in  .Monterey  Bay  along  a  13.3  km  path  from 
Pt.  Pinos  to  Marina  (Figure  2).  The  R/V  ACANIA  was  located  In  the  region 
Indicated  by  the  square.  The  aircraft  made  constant  altitude  passes  along 
the  optical  path.  Later  in  the  experiment,  several  passes  were  made  per¬ 
pendicular  to  the  path.  The  optical  comparison  was  done  on  turbulence 
(Cpj2)  and  extinction  (a  ).  The  optical  and  ship  comparison  has 
already  been  reported  (Davidson  et.  al.,  1980)  so  it  will  not  be  discussed 
here. 


1.  Aircraft  Evaluation 

The  basic  aircraft  optical  path  micrometeoro logical  measurements 
and  bulk  calculations  of  scaling  parameters  are  given  in  Table  III.  A 
more  detailed  printout  is  given  in  Appendix  A.  The  bulk  scaling  predic¬ 
tions  of  (Eq.  3  and  Fig.  1)  are  compared  with  the  turbulence  mea¬ 
surements  in  Fig.  3.  These  results  are  similar  to  those  obtained  from  the 
ship  measurements  (Davidson  et.  al.,  1980). 

2.  Aircraft  Extinction  Comparison 

A  summary  of  the  aircraft  optical  path  extinctions  is  given  in 
Appendix  B.  In  Fig.  4  the  aircraft  aerosol  comparison  with  the  optical 
measurements  is  shown.  Out  of  nineteen  comparison  runs  (three  wave¬ 
lengths  each)  only  two  disagree  by  more  than  a  factor  of  two.  For  the 
aircraft  aerosol  data,  the  average  ratio  of  extinction  for  aerosols  versus 
optics  is  1.0  +  50X,  -40%. 
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MAG  AT -80 


Comparison  of  aircraft  turbulence  Cjjj  values  (X)  and  optically  measured  Cj|j 
values  (•)  with  Cjjj  calculated  using  the  bulk  method.  The  dark  solid  line 
indicates  perfect  agreement  and  the  lighter  lines  indicate  factor  of  two  limits 
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Figure  4.  Comparison  of  optically  measured  extinction  coefficient 
and  aerosol  extinction  coefficient  from  aircraft  optical 
path  data. 
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3.  Ship  Extinction  Comparison 

The  ship  aerosol  measurements  were  made  at  anchor  along  the  opti¬ 
cal  path  or  underway  within  the  square  indicated  in  Fig.  2.  Selected  time 
series  plots  of  extinction  coefficients  are  given  in  Fig.  5a  -  5g.  Direct 
conparisons  of  size  spectral  and  optical  extinction  values  are  given  in 
Fig.  6a  -  6d. 
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Figure  5a.  Time  series  plot  of  aerosol  extinctions  coefficient 
from  optical  measurements  (lines)  and  aerosol  size 
spectra  (X  and  0). 
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Figure  5b.  Time  series  plot  of  aerosol  extinctions  coefficient 
from  optical  measurements  (lines)  and  aerosol  size 
spectra  (X  and  0).  l7 
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Figure  5c.  Time  series  plot  of  aerosol  extinctions  coefficient 
from  optical  measurements  (lines)  and  aerosol  size 
spectra  (X  and  0).  -|g 
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Figure  5d.  Time  series  plot  of  aerosol  extinctions  coefficient 
from  optical  measurements  (lines)  and  aerosol  size 
spectra  (X  and  0). 
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Figure  5e.  Time  series  plot  of  aerosol  extinctions  coefficient 
from  optical  measurements  (lines)  and  aerosol  size 
spectra  (X  and  0). 
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Figure  5f.  Time  series  plot  of  aerosol  extinctions  coefficient 
from  optical  measurements  (lines)  and  aerosol  size 
spectra  (X  and  0). 
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Figure  5g.  Time  series  plot  of  aerosol  extinctions  coefficient 
from  optical  measurements  (lines)  and  aerosol  size 
spectra  (X  and  0).  2 
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Figure  6a.  Comparison  of  aerosol  extinction  coefficient 

from  optical  and  aerosol  size  spectra  measurements 
at  \  *  0.63  u  . 
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Figure  6b.  Comparison  of  aerosol  extinction  coefficient 

from  optical  and  aerosol  size  spectra  measurements 
at  X  ■  0.84  v  . 
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Figure  6c.  Comparison  of  aerosol  extinction  coefficient 

from  optical  and  aerosol  size  spectra  measurements 
at  X  *  1 .03  u  . 
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Figure  6d,  Comparison  of  aerosol  extinction  coefficient 

from  optical  and  aerosol  size  spectra  measurements 
at  x  *  1 .06  u  . 


APPENDIX  A 


Aircraft  bulk  meteorology  and  turbulence  data  on  optical  path  In  Monterey 
Bay. 

ALT  ,  Altitude  (ft) 

PRES  ,  Pressure  (mb) 

T-ROS ,  Temperature  (cent) 

T-SIR,  Sea  surface  temperature  (cent) 

T-DEW,  Dew  point  (cent) 
q  ,  Water  vapor  mixing  ratio  (g/Kg) 

EPS  ,  Turbulence  dissipation  ratio  (m^/sec-*) 

CT2  ,  Temperature  structure  function  (K^/m^/3) 

CQ 2  ,  Water  vapor  structure  function  (ab^/m^/3) 

EXT  ,  Aerosol  extinction,  \  -0.49  (km”*) 
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APPENDIX  B 


Siamary  of  aircraft  optical  measurements  in  Monterey  Bay. 


ALT  , 

Altitude  (ft) 

PRES  , 

Pressure  (mb) 

T-ROS  , 

Tenperature  (cent) 

T-DEW  , 

Dew  point  (cent) 

CN2T  , 

Cp2  component  of  Cpj2 

:  (m-2/3) 

CN2  , 

Turbulence  value  of  %2  (m“2/3) 

E(  .63)  , 

Aerosol  extinction 

(Km-1)  at  A 

■  0.63 

E( .84)  , 

Aerosol  extinction 

(Km-1)  at  A 

-  0.84 

E(1.06) , 

Aerosol  extinction 

(Kmr^)  at  A 

-  1.06 

30 
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